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ABSTRACT 

We present a SINFONI integral field kinematical study of 33 galaxies at z ~ 3 from the AMAZE and LSD projects which are aimed 
at studying metallicity and dynamics of high-redshift galaxies. The number of galaxies analyzed in this paper constitutes a significant 
improvement compared to existing data in the literature and this is the first time that a dynamical analysis is obtained for a relatively 
large sample of galaxies at z ~ 3. 11 galaxies show ordered rotational motions (~ 30% of the sample), in these cases we estimate 
dynamical masses by modeling the gas kinematics with rotating disks and exponential mass distributions. We find dynamical masses 
in the range 2 x 10' Mq - 2 x 10"Mq with a mean value of ~ 2 x IO'^Mq. By comparing observed gas velocity dispersion with 
that expected from models, we find that most rotating objects are dynamically "hot", with intrinsic velocity dispersions of the order 
of ~ 90 km s"'. The median value of the ratio between the maximum disk rotational velocity and the intrinsic velocity dispersion 
for the rotating objects is 1.6, much lower than observed in local galaxies value (~ 10) and slightly lower than the z ~ 2 value 
(2 - 4). Finally we use the maximum rotational velocity from our modeling to build a baryonic Tully-Fisher relation at z~3. Our 
measurements indicate that z ~ 3 galaxies have lower stellar masses (by a factor of ten on average) compared to local galaxies with 
the same dynamical mass. However, the large observed scatter suggests that the Tully-Fisher relation is not yet "in place" at these 
early cosmic ages, possibly due to the young age of galaxies. A smaller dispersion of the Tuly-Fisher relation is obtained by taking 
into account the velocity dispersion with the use of the 5o5 indicator, suggesting that turbulent motions might have an important 
dynamical role. 

Key words. Galaxies: evolution - Galaxies: formation - Galaxies: fundamental parameters - Galaxies: high-redshift - Galaxies: 
kinematics and dynamics 

Cq 1. Introduction mergers are believed to play an important role for galaxy forma- 

^, , . , . ^ , . , ^ , , , tion and evolution. However, the observational evidence for the 

The dynamical properties of galaxies play a fundamental role ^^^^^^^^^ ^f rotating disks with high star formation rates at z ~ 2 
in the context of galaxy formation and evolution. The observed ^^^^ ^^^^^^ ^^^^^^^^^^ ^^ .^^^^ -^ ^j^^ ^^ ■ 

dynamics of galaxies represent a test to theoretical models and ^^^^ mechanism that drives star formation and mass assembly 

are the most direct way to probe the content of dark matter In ^^ ^.^ ^^^^^-^^ ^^-^^^ ^^ ^j 2009; Wright et al. 2009; Forster 

particular, the dynamical mass of a galaxy mferred from its rota- Schreiber et al. 2009; Cresci et al. 2009). 
tion curve is the most direct way to constrain mass and angular 
momentum of dark matter haloes, which can be compared with 

the predictions of cosmological cold dark matter models of hier- j„ recent years, our observational knowledge in this field has 

archical sti-ucture formation. In such models (Blumenthal et al. increased enormously. Many dynamical studies have been per- 

1984; Davis et al. 1985; Sprmgel et al. 2006, Mo et al. 1998) foj^ied on extended samples of z ~ 1 - 2.5 objects (Genzel et al. 



2006; Genzel et al. 2008; Forster Schreiber et al. 2006; Forster 
Schreiber et al. 2009; Cresci et al. 2009; Erb et al. 2006; Law 
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a contract with NASA. 2007; Nesvadba et al. 2008; Jones et al. 2010; Law et al. 2009; 
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Lemoine-Busserolle et al. 2009; Swinbanket al. 2007; Swinbank 
et al. 2009). 

The redshift range z ~ 3 -4 is particularly important to study 
since it is before the peak of the cosmic star formation rate (see, 
for example, Dickinson et al. 2003; Rudnick et al. 2006; Hopkins 
& Beacom 2006; Mannucci et al. 2007), when only a small frac- 
tion (~ 15%, Pozzetti et al. 2007) of the stellar mass in present- 
day galaxies has been assembled. It is also the redshift range 
when the most massive early-type galaxies are expected to form 
(see, for example, Saracco et al. 2003). Moreover, the number of 
galaxy mergers is much larger than at later times (Conselice et al. 
2007; Stewart et al. 2008). As a consequence the predictions of 
different models tend to diverge significantly at z > 2.5. 

This paper is based on two projects focused on studying 
metallicity and dynamics of high-redshift galaxies: AMAZE 
(Assessing the Mass-Abundance redshift Evolution) (Maiolino 
et al. 2008a, Maiolino et al. 2008b) and LSD (Lyman- break 
galaxies Stellar populations and Dynamics) (Mannucci & 
Maiolino 2008, Mannucci et al. 2009). Both projects use integral 
field spectroscopy of samples of z ~ 3 - 4 galaxies in order to 
derive their chemical and dynamical properties. In both projects 
we make use of data obtained with the Spectrograph for Integral 
Field Observations in the Near Infrared (SINFONl) at the Very 
Large Telescope (VLT) of the European Southern Observatory 
(ESO). Integral field spectroscopy has proven to be a powerful 
tool to study galaxy dynamics as it provides two dimensional 
velocity maps, without the restrictions of longslit studies, also 
plagued by unavoidable light losses. 

The AMAZE sample consists of ~ 30 Lyman Break Galaxies 
in the redshift range 3 < z < 4.8 (most of which at z~3.3), 
with deep Spitzer/IRAC photometry (3.6 - 8//m), an important 
piece of information to derive reliable stellar masses. These 
galaxies were observed with SINFONl in seeing-limited mode. 
The LSD sample is a representative, albeit small, sample of 10 
LBGs at z~3 with available Spitzer and HST imaging from the 
Steidel et al. (2003) catalogue. For LSD SINFONl observations 
were performed with the aid of adaptive optics in order to im- 
prove spatial resolution since this project was aimed to obtain 
spatially-resolved spectra for measuring kinematics and gradi- 
ents in emission lines. The global sample is listed in Tab.l. 

Lyman Break Galaxies are objects selected based on the Ly- 
break technique and based on their UV rest-frame blue color 
As a consequence these samples are biased both against dust 
reddeneded systems and against aged stellar populations, which 
are characterized by redder colors. Therefore, the conclusions 
inferred from our results apply only to a sub-population (about 
half) of galaxies at z~3 (Reddy et al. 2005, van Dokkum et al. 
2006, Grazian et al. 2007). 

In this paper we focus on the dynamics of AMAZE and LSD 
galaxies at z ~ 3. In Sect. 2 we present observations, data re- 
duction and the method followed to extract kinematical maps 
from integral field spectra. In Sect. 3 we present the results and 
a simple method to individuate the objects with a velocity map 
consistent with a rotating disk. In Sect. 4 we describe the "rotat- 
ing disk" model adopted in this paper, explain our "fitting strat- 
egy" to constrain model parameters and estimate their errors. In 
Sect. 4.3 we discuss the results of model fitting on the selected 
subsample of objects consistent with rotating disks. In Sect. 5 we 
discuss fit results and present the Tully-Fisher relation at z ~ 3 
(Sect. 5.2). 

In the following, we adopt a KCDM cosmology with Hq - 
10kms-^Mpc-\ 0,„ = 0.3 and Oa = 0.7. 



2. Observations and data reduction 

Complete descriptions of the AMAZE and LSD programs, of 
their observations and data reduction are presented in Maiolino 
et al. (2008a) and Mannucci et al. (2009). Here we report a brief 
summary on observations and data reduction. 

The near-IR spectroscopic observations were obtained by 
means of SINFONl, the integral field spectrometer at VLT 
(Eisenhauer et al. 2003). For AMAZE galaxies, SINFONl was 
used in its seeing-limited mode, with the 0.125" x 0.25" pixel 
scale and the H+K grism, yielding a spectral resolution R ~ 
1500 over the spectral range 1.45 - 2.41//m. For the LSD galax- 
ies, SINFONl was used with the Adaptive Optics module using 
a bright star close to the galaxy to guide the wavefront correc- 
tion system. The pixel scale used is of 0.05" x 0.10" (for all 
object except SSA22a-C6 and SSA22a-M4 for which was used 
the 0.125" X 0.25" pixel scale). The (K-band) seeing during the 
observations was generally about 0.6" -0.7". In the AO assisted 
observations the spatial resolution obtained is ~ 0.2". 

Data were reduced by using the ESO-SINFONI pipeline 
(version 3.6.1). The pipeline subtracts the sky from the tempo- 
rally contiguous frames, flat-fields the images, spectrally cal- 
ibrates each individual slice and then reconstructs the cube. 
Individual cubes were aligned in the spatial direction using the 
offsets of the position of the [OIII] or Ha line emission peak. 
Within the pipeline the pixels are resampled to a symmetric an- 
gular size of 0.125" x 0.125" or 0.05" x 0.05". The atmospheric 
absorption and instrumental response were taken into account 
and corrected by dividing with a suitable standard star. 

In table 1 we summarize all the relevant properties and ob- 
servation setups of the object presented in this paper 

In the AMAZE sample there is a subsample of 4 lensed ob- 
jects. We decided not to analyze these objects in this paper be- 
cause of the uncertainties introduced by the lensing model and 
the associated de-projection. They will be discussed in a forth- 
coming paper. 

2.1. Extraction of the gas kinematics 

We extract the kinematics of the gas by fitting the emission line 
spectrum for each spatial pixel of the cube corresponding to a 
given position on the sky. In order to improve the signal-to-noise 
ratio (hereafter S/N) we first perform a smoothing of the cube 
spatial planes, by using a gaussian filter with FWHM of 3 pixels. 
This is smaller than the spatial resolution of the observations 
which usually corresponds to a Point Spread Function with a 
FWHM of at least 4 pixels. 

In all cases, but one, we fit the profile and shift 
of [OIII]i^5007,4959A doublet. The two [OIII] lines are 
parametrized using single gaussian functions with the follow- 
ing, physically motivated constraints: the two lines are forced 
to have the same average velocity v and velocity dispersion cr, 
while their flux ratio is fixed, /^4959//^5007 = 0.33, since both the 
lines are emitted from the same upper energy level of the O^^ 
ion. For one object at z ~ 2.6 we fit the Ha line with a sin- 
gle gaussian function. Since we perform the fit of the spectrum 
within a spectral window of only ~ 0.2/iOT around the line, the 
continuum level is simply set to the mean continuum flux in this 
band. 

To automatically exclude noise fluctuations or bad pixels and 
to take into account the eff'ect of instrumental broadening we 
constrain the velocity dispersion of the lines to be larger than the 
instrumental resolution estimated from the sky emission lines. 
We define the signal-to-noise ratio (S/N) of each line as the peak 
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Table 1. Galaxy properties and observation setups. 



Object 


sample 


R.A. (" 


Dec. <" 


z/scale(kpc/")*^' 


Texp(min) 


pixel scale(") 


line used <'> 


SSA22a-M38 


AMAZE 


22:17:17.7 


+00:19:00.7 


3.294/7.48 


400 


0.125x0.25 


[OIII] 


SSA22a-C16 


AMAZE 


22:17:32.0 


+00:13:16.1 


3.068/7.65 


350 


0.125x0.25 


[OIII] 


CDFS-2528 


AMAZE 


03:32:45.5 


-27:53:33.3 


3.688/7.18 


350 


0.125x0.25 


[OIII] 


SSA22a-D17 


AMAZE 


22:17:18.9 


+00:18:16.8 


3.087/7.64 


250 


0.125x0.25 


[OIII] 


CDFA-C9 


AMAZE 


00:53:13.7 


+ 12:32:11.1 


3.212/7.54 


250 


0.125x0.25 


[OIII] 


CDFS-9313 


AMAZE 


03:32:17.2 


-27:47:54.4 


3.654/7.21 


250 


0.125x0.25 


[OIII] 


CDFS-9340<^> 


AMAZE 


03:32:17.2 


-27:47:53.4 


3.658/7.20 


250 


0.125x0.25 


[OIII] 


CDFS-11991 


AMAZE 


03:32:42.4 


-27:45:51.6 


3.661/7.25 


450 


0.125x0.25 


[OIII] 


3C324-C3 


AMAZE 


15:49:47.1 


+21:27:05.0 


3.289/7.49 


150 


0.125x0.25 


[OIII] 


DFS2237b-D29 


AMAZE 


22:39:32.7 


+ 11:55:51.7 


3.370/7.43 


250 


0.125x0.25 


[OIII] 


CDFS-5161 


AMAZE 


03:32:22.6 


-27:51:18.0 


3.66017.20 


300 


0.125x0.25 


[OIII] 


DFS2237b-C21 


AMAZE 


22:39:29.0 


+ 11:50:58.0 


3.403/7.40 


200 


0.125x0.25 


[OIII] 


SSA22a-aug96M16 


AMAZE 


22:17:30.9 


+00:13:10.7 


3.292/7.48 


250 


0.125x0.25 


[OIII] 


Q1422-D88 


AMAZE 


14:24:37.9 


+23:00:22.3 


3.752/7.13 


250 


0.125x0.25 


[OIII] 


CDFS-6664 


AMAZE 


03:32:33.3 


-27:50:7.4 


3.797/7.11 


500 


0.125x0.25 


[OIII] 


SSA22a-C36 


AMAZE 


22:17:46.1 


+00:16:43.0 


3.063/7.65 


100 


0.125x0.25 


[OIII] 


CDFS-4414 


AMAZE 


03:32:23.2 


-27:51:57.9 


3.471/7.34 


350 


0.125x0.25 


[OIII] 


CDFS-4417<4> 


AMAZE 


03:32:23.3 


-27:51:56.8 


3.473/7.34 


350 


0.125x0.25 


[OIII] 


CDFS- 12631 


AMAZE 


03:32:18.1 


-27:45:19.0 


3.709/7.17 


250 


0.125x0.25 


[OIII] 


CDFS- 13497 


AMAZE 


03:32:36.3 


-27:44:34.6 


3.413/7.39 


150 


0.125x0.25 


[OIII] 


CDFS- 14411 


AMAZE 


03:32:20.9 


-27:43:46.3 


3.599/7.25 


200 


0.125x0.25 


[OIII] 


CDFS- 16272 


AMAZE 


03:32:17.1 


-27:42:17.8 


3.619/7.24 


350 


0.125x0.25 


[OIII] 


CDFS- 16767 


AMAZE 


03:32:35.9 


-27:41:49.9 


3.624/7.24 


300 


0.125x0.25 


[OIII] 


Cosmic eye*^' 


AMAZE 


21:35:12.7 


-01:01:42.9 


3.075/<^> 


200 


0.125x0.25 


[OIII] 


Abelll689-1<^> 


AMAZE 


13:11:30.0 


-01:19:15.3 


3.770/<^> 


300 


0.125x0.25 


[OIII] 


Abelll689-2(5> 


AMAZE 


13:11:25.5 


-01:20:51.9 


4.868/<5> 


400 


0.125x0.25 


[Oil] 


Abelll689-4<5> 


AMAZE 


13:11:26.5 


-01:19:56.8 


3.038/<'> 


240 


0.125x0.25 


[OIII] 


SSA22b-C5 


LSD 


22:17:47.1 


+00:04:25.7 


3.117/7.61 


240 


0.05 X 0.1 (AO) 


[OIII] 


SSA22a-C6 


LSD 


22:17:40.9 


+00:11:26.0 


3.097/7.63 


280 


0.125 X 0.25 (AO) 


[OIII] 


SSA22a-M4<^' 


LSD 


22:17:40.9 


+00:11:27.9 


3.098/7.63 


280 


0.125 X 0.25 (AO) 


[OIII] 


SSA22a-C30 


LSD 


22:17:19.3 


+00:15:44.7 


3.104/7.62 


240 


0.05 X 0.1 (AO) 


[OIII] 


Q0302-C131 


LSD 


03:04:35.0 


-00:11:18.3 


3.240/7.51 


240 


0.05 X 0.1 (AO) 


[OIII] 


Q0302-C171 


LSD 


03:04:44.3 


-00:08:23.2 


3.34/7.44 


240 


0.05 X 0.1 (AO) 


[OIII] 


DSF2237b-D28 


LSD 


22:39:20.2 


+ 11:55:11.3 


2.938/7.75 


240 


0.05 X 0.1 (AO) 


[OIII] 


Q0302-M80 


LSD 


03:04:45.7 


-00:13:40.6 


3.416/7.39 


240 


0.05 X 0.1 (AO) 


[OIII] 


DSF2237b-MD19 


LSD 


22:39:21.1 


+ 11:48:27.7 


2.616/7.99 


200 


0.05 X 0.1 (AO) 


Ho- 



(1) 

(2) 
(3) 
(4) 
(5) 



12000. 

Redshifts are measured from the observed [OIII] line wavelength. 

[OIII] denotes the [OIII]^^ 5007,4959A doublet; Ho- the Ha A 6563 A line. 

The object is in the same field of view of the object on the previous line 

Lensed objects not analyzed in this paper. The spatial scale (pc/") depends on the lensing model. 



of the line model divided by the r.m.s. of the spectrum estimated 
in regions with no emission lines. Finally, we reject all fits where 
the S/N is lower than 2. 

Uncertainties on measured fluxes, velocities and velocity 
dispersions are obtained from the formal errors on best fit pa- 
rameters. Such errors are computed after two iterations of the 
line fitting procedure: in the second fit we adopt pixel by pixel 
errors equal to the r.m.s. of the residuals in the first fit. 



3. Results 

In Fig. 1 to 6 we present the kinematical maps for all the objects 
described in tab. 1 . This paper is dedicated only to the study of 
galaxy dynamical properties while other studies (e.g. integrated 
properties, metallicity, etc.) have been or will be presented in 
other dedicated papers (Maiolino et al. 2008a, Mannucci et al. 
2009, Mannucci et al. 2010, Cresci et al. 2010 and Troncoso et 
al. 2010, in prep.). 



In the majority of the AMAZE sample (17 out of 23 objects, 
i.e. ~ 72%, excluding the four lensed objects) the morphology 
of line flux maps is simple, characterized by single peak distri- 
butions with typical dimensions of ~ 1-1.5" and relatively low 
ellipticities of the isophotes. These morphologies are likely the 
consequence of beam smearing since the typical seeing in this 
data is ~ 0.6 - 0.7" FWHM (see below for a more accurate esti- 
mation). However, in some cases we find a secondary clump or 
asymmetric extended emission (e.g. in SSA22-M38). 

In the galaxies from the LSD sample we find more com- 
plex morphologies which are likely the consequence of the AO- 
assisted observations providing higher spatial resolution (~ 0.2" 
FWHM) which allows us to spatially resolve more complex 
structures, but lower sensitivity to extended sources given the 
higher surface brightness detection threshold. 

There are three cases of close pairs of objects with simi- 
lar brightnesses. CDFS-9313 and CDFS-9340 have a projected 
plane of the sky separation of ~ 1 .0" (~ Ikpc at the average red- 
shift of the sources); CDFS-9340 is redshifted by ~ 2EQkms^^ 
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Fig.l. Kinematical maps for objects SSA22A-M38, SSA22A-C16, CDFS-2528 and SSA22A-D17 (from the top respectively). Respectively, from 
the left: flux, velocity and sigma map. the X - Y coordinates are in arcseconds referred to an arbitrary object centre position. The North direction 
is the positive Y axis. The vertical color bars are in arbitrary units for the flux map and in km .s"' for the velocity and sigma maps. Overplotted on 
the flux map the continuum flux distribution (brown isophotes) for the object whenever we detect a continuum flux component with sufficient S/N. 



with respect to CDFS-9313, the brightest of the pair CDFS- 
4414 and CDFS-4417 have a projected separation of ~ 1.2" 
(~ 9kpc); CDFS-4414 is blueshifted by ~ llQkms-^ with re- 
spect to the brighter companion. SSA22a-C6 and SSA22a-M4, 
from the LSD sample, have a projected separation of ~ 1.6" 
(~ Ukpc) and SSA22a-M4 is redshifted by ~ 90/tm.r' with re- 
spect to the brighter companion. 

Inspection of the velocity maps in Fig. 1 to 6 reveals a num- 
ber of objects with a clear and regular velocity gradient (12 out 



of 22 objects in the AMAZE sample and only one out of 9 in the 
LSD sample). The presence of such velocity gradients is an in- 
dication of possible rotating-disk kinematics. The identification 
of rotating objects and their distinction from complex kinemat- 
ics cases or mergers is important for the comparison with the- 
oretical models. Some authours (Forster Schreiber et al. 2009, 
Cresci et al. 2009) used the technique developed by Shapiro et al. 
(2008) based on kinemetry (Kxajnovic et al. 2006) to quantify 
asymmetries in both the velocity and sigma maps in order to 
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Fig. 2. Kinematic maps for the objects (respectively from the top) CDFa-C9, CDFS-9313 (in the same field of view the fainter source at the 
North- West is CDFS-9340), CDFS-11991, 3C324-C3 and DFS2237b-D29. Panels, axis and color bar indications same as in Fig. 1. 
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Fig. 3. Kinematic maps for the objects (respectively from tlie top) CDFS-5161, DFS2237b-C21, 
Panels, axis and color bar indications same as in Fig. 1. 
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Fig. 4. Kinematic maps for the objects (respectively from the top) SSA22A-C36, CDFS-4414/4417(CDFS-4417 is the brighter object at the north), 
CDFS- 1263 1,CDFS- 13497 and CDFS-14411. Panels, axis and color bar indications same as in Fig. 1. 
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Fig. 5. Kinematic maps for the objects (respectively from the top) CDFS- 16272, CDFS-16767, SSA22b-C5, SSA22a-C6/M4 (SSA22a-M4 is the 
fainter northern object) and SSA22a-C30. Panels, axis and color bar indications same as in Fig. 1. 
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Fig. 6. Kinematic maps for tiie objects (respectively from the top) Q0302-C131, Q0302-C171, DSF2237b-D28, Q0302-M80 and DSF2237b- 
MD19. Panels, axis and color bar indications same as in Fig. 1. 
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empirically differentiate between rotating and non rotating sys- 
tems. This technique is not useful when applied to the data pre- 
sented in this paper due to the low S/N that does not allow us to 
recover with sufficient accuracy the needed kinemetry parame- 
ters. We assess the presence of such velocity gradients by using 
a different approach. We fit the velocity map with a plane (i.e. 
V - ax + by + vq) and classify the object as "rotating" or "not 
rotating" if the velocity map is well fitted (or not) by a plane. In 
detail we use the x^ statistic to accept or reject the hypothesis 
that the velocity map is well fitted by a plane: we accept the hy- 
pothesis if there is a probability lower than the 4.6% (adopting a 
X^ distribution with the correct number of degree of freedom) to 
obtain ax^ value greater than the value obtained. 

For objects with velocity maps well fitted by a plane we cal- 
culate its inclination relative to the X- F plane (i.e the magnitude 
of the velocity gradient). We confirm its "rotating" classification 
only if this inclination is inconsistent with zero at least a signifi- 
cance level of Act. On the other hand a plane inclination consis- 
tent with zero might correspond either to genuinely non rotating 
object or to a rotating one seen face-on. In this case we cannot 
identify it as a rotating disk. In conclusion we classify the ob- 
served galaxies with the following criteria. 

- If the velocity map shows a non zero gradient from plane 
fitting, we classify the object as "rotating". 

- If the velocity map is well fitted by a plane but its inclination 
is consistent with zero within 4cr, we label the object as "not 
classifiable". 

- If the velocity map cannot be fitted with a plane, we classify 
the object as "not rotating". 

We choose to fit a "plane" to the velocity gradient because, 
due to the limited S/N, we are sampling preferentially the linear 
part of the rotation curve of our targets. By using a plane we can 
adopt a simple model with few free parameters and therefore we 
can rely on more simple statistics. 

In principle we could have used also velocity dispersion 
maps to classify galaxies as "rotating" or "non rotating": in fact, 
a circularly symmetric dispersion map is the expected signature 
of spatially unresolved rotation at the center of a rotating disk. 
However, we decided not to use them principally because we 
would have missed galaxies with regular rotation patterns in ve- 
locity, but whose velocity dispersion maps are distorted by turbu- 
lent and/or non-gravitational motions, as verified "a posteriori" 
with our modeling. 

In Fig. 7 we show an example of our method. We plot the 
projection of the velocity map along the direction (s) of max- 
imum inclination of the fitted plane (i.e. the fitted plane seen 
from its edge). We also show the 1,2,3, and 4cr slopes of the 
plane inclination. 

The classification for all the objects is reported in table 2. 

Note that the object SSA22a-M38 is classified as "rotating" 
even if its velocity map is not fitted by a plane (see the x^ value 
in tab. 2). In this case the smooth velocity gradient is very clear 
and the S/N is high; the fit with a simple plane fails because of 
the flattening of the velocity field at the "red" and "blue" edges. 

The next step is to model the gas kinematics of "rotating" 
galaxies to relate observed kinematics with the mass distribution. 

4. Kinematical modeling 

4. 1 . The model 

The adopted dynamical model assumes that the ionized gas is 
circularly rotating in a thin disk. We neglect all hydrodynamical 




-2 2 

3 (pixel) 

Fig. 7. Example of the method (applied to object CDFA-C9) for assess- 
ing the presence of a smooth gradient in the velocity map. The solid line 
shows the fitted plane seen from the edge, while the filled points show 
the observed velocity map on the projection. The dashed lines represent 
the 1,2,3, and 4cr on the plane inclination. 

Table 2. Objects observed. 



AMAZE Objects 


X^ld.o.f. 


Classification 


SSA22a-M38 


3057/176 


Rotating 


SSA22a-C16 


182/136 


Rotating 


CDFS-2528 


61.6/52 


Rotating 


SSA22a-D17 


20.6/42 


Rotating 


CDFA-C9 


157/142 


Rotating 


CDFS-9313 


63.8/91 


Rotating 


CDFS-9340 


42.5/29 


Rotating 


CDFS- 11991 


226/99 


Not rot. 


3C324-C3 


45.3/76 


Rotating 


DFS2237b-D29 


229/107 


Not rot. 


CDFS-5161 


2.00/12 


Not clas. 


DFS2237b-C21 


725/141 


Not rot. 


SSA22a-aug96M16 


680/119 


Not rot. 


Q1422-D88 


478/127 


Not rot. 


CDFS-6664 


15.9/40 


Not clas. 


SSA22a-C36 


102/86 


Not clas. 


CDFS-4414 


259/91 


Not rot. 


CDFS-4417 


117/94 


Not clas. 


CDFS-12631 


135/70 


Not rot. 


CDFS- 13497 


3133/91 


Not rot. 


CDFS- 144 11 


62.7/68 


Rotating 


CDFS-16272 


48.4/82 


Not clas. 


CDFS-16767 


63.0/46 


Rotating 


LSD Objects 


SSA22b-C5 


216/138 


Not rot. 


SSA22a-C6 


101/54 


Not rot. 


SSA22a-M4 


21.8/36 


Not clas. 


SSA22a-C30 


19334/243 


Not rot. 


Q0302-C131 


62.2/95 


Rotating 


Q0302-C171 


170/102 


Not clas. 


DSF2237b-D28 


1848/203 


Not rot. 


Q0302-M80 


206/87 


Not rot. 


DSF2237b-MD19 


268/95 


Not rot. 
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effects, therefore the disk motion is entirely determined by the 
gravitational potential. 

In principle turbulent motions, often observed in high z ro- 
tating disk galaxies as well as in our own data, could provide 
significant dynamical support. Indeed Epinat et al. (2009) ap- 
ply an asymmetric drift correction assuming that turbulent mo- 
tions support mass. However, while stellar chaotic motions al- 
ways support mass because they have gravitational origin, this 
is not always true for gas motions. The gas velocity dispersion 
might be strongly increased because of non gravitational mo- 
tions (like, e.g., in winds driven by starburst activity), and in 
that case it would not be directly linked to the dynamical mass. 
Therefore we adopted a conservative approach and neglected the 
contribution of the gas velocity dispersion. In any case we will, 
at most, underestimate the dynamical mass by a factor lower than 
~ 2 at the disk scale length (Epinat et al. 2009). 

The galaxy gravitational potential is generated by an expo- 
nential disk mass distribution with surface density given by: 



2(r) = Eo e 



-rjr. 



(1) 



The rotation curve of a thin disk with due to such mass distribu- 
tion is (Binney & Tremaine 1987): 



Ycirf = \TiGY.^r,y\h{y)KQ{y) - h{y)K^{y)\ 



(2) 



where y - rjlr^ and /„, K„ are the modified Bessel functions of 
the first and second kind. 

By defining the "dynamical mass" of the galaxy as the total 
mass enclosed in a 10 kpc radius (e.g. Cresci et al. 2009), we can 
finally write 



Mdyn 



10"'Mfr 







yc(r) = 6.56 km s 

where 

HyWoiy) - h{y)Kx{y) 



(^)->^' 



A{y) 



/o(3'io)^o(}'io) - hiyioWiiyia) 



0.5 



(3) 



(4) 



and yio corresponds to y computed for r = 10 kpc. 

The velocity along the line of sight for a given position on 
the sky is derived from Vdr) by taking into account geometrical 
projection effects. 

Another important element of the model is the intrinsic flux 
distribution (hereafter IFD) of the emission line, since it acts as 
a weighting function in the calculation of the observed velocities 
and velocity dispersions when taking into account instrumental 
effects. Following Cresci et al. (2009) the IFD in the disk plane 
is modeled with an exponential function. 



/(r) = he 



-r/ro 



(5) 



where r is the distance from the disk center and fq is the scale 
radius. Such IFD is then projected onto the plane of the sky. 
Several objects in our sample have surface brightness distribu- 
tions clearly different from those of exponential disks. However, 
given the low S/N and poor spatial resolution of our data, we de- 
cided to keep the IFD modeling as simple as possible (see also 
Cresci et al. 2009). 

In computing the observed velocity shift and velocity disper- 
sion we take into account the instrumental beam smearing. The 
spatial PSF is modeled with a two-dimensional gaussian func- 
tion with full width half maximum FWHMpsF, sampled over 
the detector x, y pixels. The instrumental spectral response is 
also described by a Gaussian function with sigma cr,„5,. Finally, 



we calculate the model values for the three observed quantities: 
Fm^dei(xi,yj), Vmodei{xi,y j) and o-,nodei{xi,y j) where Xi,yj repre- 
sent a spatial pixel in the data cube. For a detailed description of 
the model refer to appendix B of Marconi et al. (2003). 
Summarizing, the model parameters are: 



Xc,yc 



I 

Mdyn 
re 



coordinates in the plane of the sky of the disk dy- 
namical centre. 

position angle (PA) of the disk line of nodes, 
inclination of the disk, 
dynamical mass. 

characteristic radius of the exponential disk. 
systemic velocity of the galaxy. 



By using these parameters values we derive Vmax that is the 
maximum velocity of the rotation curve. We will use this param- 
eter in Sect. 5.2 to build the Tully-Fisher relation for our data 
sample. 

4.2. Fit Strategy 

4.2.1. Preliminary steps 

The first preliminary step of the model fitting is to estimate the 
spectral and spatial resolution of the observations, character- 
ized by the corresponding cTspec and FWHMpsF of the gaussian 
broadening functions. 

The spectral resolution {(Tspec) is simply estimated by using 
the profile of isolated lines in sky spectra. 

Estimating the spatial resolution {FWHMpsp) is more com- 
plex since we miss datacubes of stars obtained with the same 
instrumental setting at the same time and at a similar position 
on the sky as the object datacube. We have therefore devised a 
method to estimate FWHMpsp directly from the datacubes of 
the objects showing rotating-disk kinematics. 

In the assumption that the kinematics is well approximated 
by a rotating disk we expect an intrinsic line of sight velocity 
field on the sky plane as shown in Fig. 8, which illustrates the 
iso-velocity contours projected on the sky of a simulated galaxy 
with mass 1O'"M0 (the so called "spider diagram"). The gas 
with line of sight velocity in a given velocity bin lies in the locus 
delimited by two subsequent isovelocity contours. We can ob- 
serve that the spatial extent of such region is always lower along 
the line of nodes direction (the X-axis) than in the perpendicular 
direction, and the central velocity bins show the minimal extent. 
For this reason we expect that an image of the object in the cen- 
tral velocity bins will be spatially unresolved along the line of 
nodes direction (see the example of Fig. 8: the central velocity 
bin of 30km ,? ' has an extent along the line of nodes direction 
lower than ~ 0.1", well below the typical value of our spatial 
resolution of ~ 0.6" FWHM). 

In practice, we select the spatial planes of the object datacube 
with wavelength close the line centroid, and for each of these 
spatial planes we fit the line surface brightness with a two di- 
mensional gaussian function. According to what is stated above, 
we take the minimum FWHM value from these two dimensional 
gaussians as an estimate of FWHMpsp- 

We could test our method with one object, CDFa-C9, for 
which there is a star in the SINFONI field of view. The FWHM 
of this star is consistent within ~ 0.1" with our estimate of the 
spatial resolution. We also verified that our estimates are in good 
agreement with the FWHM of the telluric standard stars ob- 
served after the science exposures, although with different air- 
masses. Moreover, in the following, we will test how systematic 
errors on FWHMpsp can affect the best fit model parameters. 
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log(Mdyn/Msun)=10 Scale=7482pc/arcsec i=60deg 
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Fig. 8. Velocity field on the sky plane ("spider diagram") for a simulated 
rotating disk at redshift z ~ 3.3 with mass of IO'^Mq, an exponential 
mass profile with characteristic radius of ~ l.lkpc, inclination of 60° 
and velocities binned in steps of 30fcm i-"'. 



The second preliminary step is the determination of the in- 
trinsic flux distribution (IFD) to be used in the computation of 
the velocity maps. The exponential IFD in the disk plane is pro- 
jected onto the plane of the sky, convolved by the PSF, averaged 
over the detector binsize and matched to the observed flux map 
by minimizing the following ;^f^ 



Xflux 



z 

L Ui 



F{xi,yj)-F,nodei{xi,yj) 



^F(xi,yj) 



(6) 



the free parameters in this fit are the center position of the 
IDF in the plane of the sky {xc,yc), the scale radius ro and the 
position angle and inclination of the elliptical isophotes (/q, the 
flux scale, is only a normalization factor). 

From the fit of the IFD we also estimate some parameters 
needed in the model. As explained in the following section, 
we set the disk center position equal to the IFD center position 
(xc,yc) and the mass distribution scale radius r^ equal to the IFD 
scale radius ro. The inclination and position angle of the disk are 
instead determined from the fit of the velocity maps, but we use 
the estimates from the IFD fit as first guesses. 

4.2.2. Fitting of velocity maps 

To estimate the free model parameters we minimize the follow- 
ing ^2. 



2 



z 



v(xi,yj) - v„,odei(xi,yj) 



n2 



Av{xi,yj) 



(7) 



a different approach in which some of these parameters are con- 
strained a priori. 

The mass distribution should be traced by the continuum flux 
distribution, therefore the disk center position should be identi- 
fied by the continuum peak. Only in some objects we can detect 
continuum emission with high enough S/N to locate the disk cen- 
ter In such cases we use that position in the fit allowing a vari- 
ation of ±0.1" consistent with measurement errors. For all other 
objects we use the position inferred from the fitting of the emis- 
sion line flux maps. For those objects where we detect a contin- 
uum component the average distance between line and contin- 
uum peaks is ~ 0.25". Therefore, when we do not detect any 
continuum, we allow for a variation of the disk center position 
of ±0.25". 

The scale radius of the mass distribution, r^ , cannot be de- 
rived from the velocity maps. Following Cresci et al. (2009), we 
fix /•(, to the value of ro, the scaling radius of the continuum flux 
map. When this is not detected, we adopt the scale radius of the 
emission line flux map. Estimating the scale radius re by fitting 
the IFD with a particular model (i.e. exponential) can of course 
lead to systematic errors. Therefore we will analyze for every 
object how systematic errors on r^ affect dynamical mass val- 
ues. 

Another potential source of systematic errors is the use of 
the line flux map to estimate disk center position and scale ra- 
dius when no continuum is detected. Therefore we also evalu- 
ated how this choice affects the dynamical mass: for the objects 
where the continuum is detected we performed the fit of the ve- 
locity maps using the disk center position and scale radii esti- 
mated either from continuum and line flux maps. We find an 
average variation of 0.2 dex for the best fit M^,,, values. This can 
be considered an estimate of the systematic error associated to 
the use of the line flux map instead of the continuum one. 

Since the disk inclination is coupled with the total dynamical 
mass, we decide to keep it as a fixed parameter in the fitting pro- 
cedure to avoid convergency problems with the;i'^ minimization. 
Its best value and confidence interval is then identified by using 
a grid of / values and performing the fit for each of them. In this 
way we construct the x^ curve as a function of / which allows us 
to identify both the best / value and the confidence intervals (for 
details see Avni 1976). In the cases where the inclination is not 
constrained by the fit, we set its value to 60°, which represents 
the mean value for a population of uniformly randomly oriented 
disks. 

4.2.3. Tine intrinsic velocity dispersion 

Our model computes the observed line velocity dispersion by 
taking into account the unresolved rotation, the intrinsic instru- 
mental broadening ((Tspec) and the broadening due to the finite 
size of the spatial pixels and of the spatial PSF. However, the 
disk can be characterized by an intrinsic velocity dispersion due 
to turbulent and non-gravitational motions in general. We esti- 
mate such intrinsic velocity dispersion of the gas as 

0",„,(X;, yjf = max[ (o-otservedixi, y jf - CTmodAXi, yjf) ', 0] (8) 

where o-ohservedixi,yj) is the observed map of velocity dispersion 
and cr,nodei(xi, yj) is that computed with the model which best fits 
the velocity map. 



In principle, all model parameters (disk center position, orien- 
tation and inclination, systemic velocity and dynamical mass) 
should be left free to vary. However, the moderate S/N of the 
data combined with the limited spatial coverage suggest to adopt 



4.2.4. Error estimates 

To estimate the errors on the best fitting model parameters we 
use the bootstrap method (Efron & Tibshirani 1994). The dataset 
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of each galaxy consists of n data points each characterized by 
spatial position, flux, velocity and velocity dispersion. We ran- 
domly extract from each dataset a subsample of n data points. 
Due to the random extraction, each subsample will have some 
data points that are replicated a few times and some data points 
that are entirely missing. We perform the fit on 100 subsamples 
and we then estimate errors on parameters by taking the stan- 
dard deviation of the best fit values which are usually normally 
distributed. 



4.3. Fit Results 

The fit procedure outlined in the previous section is applied to 
the subsample of "rotating" objects in Tab. 2. In the following 
we present and discuss the results of the fit for each of these ob- 
jects. The observed, model and residual velocity maps are pre- 
sented in the top panels of Figs. 9-19. The bottom panels of the 
same figures show the observed and model and intrinsic velocity 
dispersion maps. The best fit parameters are presented in Tab. 3. 



SSA22a-M38 (Fig. 9). In this object, as shown in Fig. 1, we de- 
tect continuum emission which we use to estimate the po- 
sition of the disk center. The inclination is well constrained 
by the observations to / - 83° ^^ at the Icr confidence level 
indicating that the disk is seen edge-on. This is a clear con- 
sequence of the fact that the observed isovelocity contours 
do not show any curvature as expected when the rotating 
disk is seen at intermediate inclinations. The best fitting dy- 
namical mass is then little affected by systematic errors due 
to uncertainties on the disk inclination. The FWHM of the 
spatial PSF estimated as described in Sec. 4.2.1 is 0.5". To 
evaluate how an incorrect estimate of this parameter can af- 
fect the Mdyn value we repeated the fit using FWHM in the 
range 0.4 - 0.6". The value of Mdyn differ in this case by 
only 0.07 dex. The characteristic radius r^ of each object de- 
pends on the value of the adopted PSF FWHM (e.g. setting 
a narrower PSF corresponds to an intrinsically broader flux 
distribution, with a larger r^,). Therefore we have also evalu- 
ated how an incorrect estimate of the scale radius r^ affects 
the Mdy„ best fit value. We repeated the fit using r^ different 
by ±30% but keeping the PSF width fixed: the value of Mdvn 
changes by only ±0.1 dex. 

By adding in quadrature all the systematic errors we obtain 
the values for the systematic errors on the dynamical masses 
reported in Tab. 3. 

In the bottom right panel of Fig. 9 we present the intrinsic 
velocity dispersion that is not accounted for by the ordered 
rotation of the gas. Most of the map is characterized by a 
very low intrinsic dispersion with average value of (cr,„,) ~ 
5Qkms^^. This means that most of the rotating disk is dy- 
namically "cold". However, a North-East region is appar- 
ently characterized by high turbulence (cr,„, ~ 25Qkms^^), 
although it is not clear whether the large velocity dispersion 
at this location is a consequence of low S/N of the emission 
lines. There is the possibility that the weak clump observable 
in the flux map at the south of the principal peak (Fig. 1) 
could be a distinct object. In our analysis we assume that 
all emission comes from an unique disk, but if the southern 
clump is really a distinct object we have to exclude its emis- 
sion from the disk. However our assumption is strengthened 
by the continuum flux map that do not show any secondary 
clump, hence the secondary sourthern clump is likely a star 
forming region within the disk. 



SSA22a-C16 (Fig. 10). In this object we detect continuum 
emission that allows us to fix the center of rotation, but the 
S/N does not allow us to estimate the scaling radius, which 
is then inferred from the flux distribution of the emission 
line. The inclination is constrained to be / = 60° ^^3^ at 
the Icr confidence level resulting in a systematic error on 
the mass of +0.45 - 0.03 dex. The PSF width estimated 
is FWHMpsF = 0.5". Considering a range of variation of 
FWHMpsF - 0.4 - 0.6", the value of Mdy„ has a systematic 
error of 0.15 dex. A variation of 30% of the disk scale radius 
produces a variation of O.ldex on the best fit M^vn value. By 
combining all the contributions to the systematic error we 
obtain a total error of +0.48 - 0.18 dex. The intrinsic ve- 
locity dispersion maps in the bottom right panel of Fig. 10 
indicates that this object is turbulent, with a quite constant 
value of (Tin, ~ lOOkms^K 

CDFS-2528 (Fig. 11). The continuum emission detected in this 
object can only be used to constrain the position of the center 
of rotation. The inclination is only constrained to be / < 70°, 
therefore we can only provide a lower limit on the dynami- 
cal mass of log (Mdyn/MQ) > 10.65. In this case, when incli- 
nation is partially or totally unconstrained, we estimate the 
mass by adopting the fiducial value of / = 60°. The PSF es- 
timated width is FWHMpsF = 0.7". Allowing a variation in 
the range FWHMpsr - 0.6 - 0.8" the value of Mj,,„ varies 
by 0.07 dex. A variation of 30% of the disk scale radius pro- 
duces a variation on M^y„ of O.ldex. The intrinsic dispersion 
is quite low on average {{(Tint) ~ 2Qkms^^) with an excess of 
~ lOOkms'^ in the East part. 

SSA22a-D17 (Fig. 12 ). For this object we do not detect any 
continuum emission. The inclination is constrained to be 
/ > 70° with a best fit value of 85°. The systematic er- 
ror on the dynamical mass due to the inclination is then 
+0.02 - 0.13 dex. The PSF width is FWHMpsf = 0.65". 
Varying FWHMpsf by ±0.1" we find a negligible system- 
atic error on the dynamical mass. Varying by 30% the disk 
scale radius produces a variation on M^yn of O.O^dex. As 
previously noted, since we do not detect any continuum in 
this object we have to consider an additional systematic er- 
ror of 0.2 dex due to the use of the line flux distribution. The 
disk is quite turbulent with an average intrinsic dispersion of 
{cint) ~ l3Qkms^^ and a southern peak of ~ 250kms^K 

CDFa-C9 (Fig. 13). For this object we detect the continuum, 
but we can only constrain the disk center position. The in- 
clination is constrained to be / > 20° with a best fit value 
of 60°. The systematic error on the dynamical mass due to 
the inclination is then +0.8 - 0.1 dex. The PSF width is 
FWHMpsF = 0.5". Varying FWHMpsf by ±0.1" we find 
a systematic error on Mj,,„ of 0.22 dex. Varying by 30% 
the disk scale radius produces a variation on Mdyn of only 
0.02dex, therefore we can neglect this contribution to the 
systematic error The disk is quite turbulent with an average 
intrinsic dispersion of (cr,„() ~ l2Qkms^^. 

CDFS-9313 (Fig. 14). (In the same field as CDFS-9340). The 
continuum is not detected. We set the inclination to / = 
60° because it is totally unconstrained. The PSF width is 
FWHMpsF = 0.7". Varying FWHMpsf by ±0.1" we find a 
Mdyn systematic error of 0.88 dex. Varying by 30% the disk 
scale radius produces a variation on Mdyn of only Q.Q2dex, so 
we can neglect this contribution to the systematic error. The 
intrinsic velocity dispersion is quite constant over the entire 
map with an average value of {o-i„,} ~ lOOkms^^. 

CDFS-9340 (Fig. 15). (In the same field as CDFS-9313). The 
continuum is not detected. The inclination best fit value 
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Fig. 9. Fit results for the object SSA22A-M38. Respectively on top panels from the left: data, model and residuals map. On bottom panels from the 
left: data, model and "intrinsic sigma" map. The plus sign on the middle panels represents the position of the model dynamical center. The X - Y 
coordinates are in arcsec referred to an arbitrary object center position. The North direction is the positive y axis. The vertical color bars are in 
km .s"'. 
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Fig. 10. Fit results for the object SSA22A-C16. Respectively on top panels from the left: data, model and residuals map. On bottom panels from 
the left: data, model and "intrinsic sigma" map. Axis, symbols and color bars indication same as in Fig. 9. 



is / - 70°, but it is unconstrained at the Icr confidence 
level. The PSF width is the same estimated for CDFS-9313. 
Varying FWHMpsF by +0.1" we find a Mdy„ systematic er- 
ror contribution of 0.1 dex. Varying by 30% the disk scale 
radius produces a variation on M,/y„ ofO.ldex. The intrinsic 
dispersion has an average value of {cr j„f} ~ 43kms^\ 
The small separation of this object from CDFS-9313 either 
in position (~ 1,0" - 1 .2kpc) and in redshift (~ 2%Qkms^^) 



suggests to consider them as an interacting pair The most 
massive object is in this case the least luminous, CDFS-9340 
(see tab. 3). We estimate a mass ratio for the pair of ~ 10, 
indicative of a minor merging event with a relatively large 
separation of the two objects (i.e. ~ l.lkpc compared to the 
disk scale radii of QJkpc and 1 Akpc), consistently with the 
non perturbed rotating disk kinematics observed in the two 
objects. 
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Fig. 11. Fit results for the object CDFS-2528. Panels, axis, symbols and color bars indication same as in Fig. 9. 



Velocity map 



Model 



Residuals 



1.0 

0.5 

0.0 

-0.5 

-1.0 





1.0 


1 


y 


0.5- 

0.0- 

-0.5 

-1.0 


> 
1 



-1.0-0.5 0.0 0.5 1.0 -1.0-0.5 0.0 0.5 1.0 

Sigma map Model 




1.0 






60 
40 


0.5 




^ 


20 


0.0 




1 J' 





-0.5 




■' 


-20 


-1.0 






-40 




-1.0 


-0 5 0-0 0-5 10 






Intrinsic Sigma 






' 


I 


250 


1.0 




1 


200 


0.5 




.-Jh 


150 


0.0 




^9 








M^ 1 


100 


-0.5 






50 


-1.0 




1 





-1.0-0.5 0.0 0.5 1.0 -1.0-0.5 0.0 0.5 1.0 -1.0-0.5 0.0 0.5 1.0 

Fig. 12. Fit results for the object SSA22a-D17. Panels, axis, symbols and color bars indication same as in Fig. 9. 



3C324-C3 (Fig. 16). In this object we do not detect any contin- 
uum emission. The inclination is constrained to be 30° < / < 
70° with a best fit value of 60° providing a systematic error 
on the dynamical mass of +0.43 - 0.05 dex. The PSF width 
is FWHMpsF = 0.7". Varying FWHMpsf by +0.1" we find 
a negligible Mjv« systematic error contribution (0.001 dex). 
Varying by 30% the disk scale radius results in a variation 
of Mdyn by Q.ldex. We note that the observed velocity field 
although showing a regular gradient is quite different from 
what expected in the case of a rotating disk. The intrinsic ve- 
locity dispersion is quite constant over the entire map with 
an average value of (o",„,) ~ l25kms^K 



CDFS-14411 (Fig. 17). For this object we do detect a contin- 
uum emission and use it to constrain the disk center posi- 
tion and scale radius. The inclination is constrained to be 
80° < / < 88° with a best fit value of 85°. The dynamic 
mass systematic error due to the inclination is then negli- 
gible (+0.04 dex). The PSF width is FWHMpsf = 0.6". 
Varying FWHMpsf by ±0.1" we find a Mdy„ systematic er- 
ror contribution of 0.17 dex. Varying by 30% the disk scale 
radius varies Mjy,, by Q.lldex. The intrinsic velocity disper- 
sion rises in the center (cr,„, ~ llQkms^^). The mean value 
over the entire map is {o-j,,,} ~ 65kms^^. 
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Fig. 13. Fit results for the object CDFa-C9. Panels, axis, symbols and color bars indication same as in Fig. 9. 
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Fig. 14. Fit results for the object CDFS-9313. Panels, axis, symbols and color bars indication same as in Fig. 9. 



CDFS-16767 (Fig. 18). In this object we detect the continuum 
emission and use it to constrain the disk center position and 
scale radius. The inclination best value is / = 15°, but it is 
unconstrained at the Icr confidence level. The PSF width is 
FWHMpsF = 0.55". Varying FWHMpsf by +0.1" we find 
a Mdy„ systematic error contribution of 0.21 dex. Varying of 
30% the disk scale radius produces a negligible variation on 
Mdy„ (0.04 dex). The disk is characterized by high turbulence 
with an average intrinsic dispersion of {o-i„,} ~ IQOkms'^. 

Q0302-C131 (Fig. 19). This is the only rotating object in the 
LSD sample. In this object we do not detect any contin- 
uum emission. The inchnation is constrained to be / < 60°, 



with a best fit value of 30° therefore we can only provide a 
lower limit on the dynamical mass of log (Mdvn/MQ) > 9.59. 
The PSF width is FWHMpsf = 0.3". Varying FWHMpsf 
by ±0.1" we find a M^,,, systematic error contribution of 
0.88 dex. Varying by 30% the disk scale radius results in a 
negligible variation on Mdyn (0.01 dex). The disk has a quite 
low average intrinsic dispersion {{cri„,} ~ 43kms^^). 
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Fig. 15. Fit results for the object CDFS-9340. Panels, axis, symbols and color bars indication same as in Fig. 9. 
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Fig. 16. Fit results for the object 3C324-C3. Panels, axis, symbols and color bars indication same as in Fig. 9. 



5. Discussion 

5.1. Rotating galaxies, turbulence and dynamical masses. 

In the AMAZE sample ~ 40% of the galaxies show a smooth ve- 
locity gradient (10 out of 24) consistent with rotating disk kine- 
matics. In the LSD sample only ~ 13% of the galaxies show 
a smooth velocity gradient (1 out of 9). The lower fraction of 
rotating objects in the LSD sample is likely due to the lower 
sensitivity of these observations to the outer, low surface bright- 
ness region (which dominate the "rotation signal") because of 
the finer pixel scale, which makes the detector read-out noise 
more important relative to the AMAZE observations performed 



with larger pixel scale. Moreover, the AMAZE sample also tend 
to have a higher fraction of bright and redder targets (generally 
translating into higher masses), which may imply a larger frac- 
tion of "settled" massive systems. We note that the fraction of ro- 
tating galaxies is fully consistent, within statistical fluctuations, 
with the fraction of rotating galaxies (~ 33%) found in the SINS 
sample at z~ 1.5-2. 5 (Forster Schreiber et al. (2009)). This is 
particularly interesting, given that the latter sample includes not 
only rest-frame UV selected galaxies, but also near- and mid-IR 
selected, and that the fraction of rotating galaxies does not ap- 
pear to depend on the selection criteria. Taken at face value, the 
comparison between the SINS analysis at z~2 and the AMAZE 
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Fig. 17. Fit results for the object CDFS-14411. Panels, axis, symbols and color bars indication same as in Fig. 9. 
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Fig. 18. Fit results for the object CDFS-16767. Panels, axis, symbols and color bars indication same as in Fig. 9. 



analysis at z~3.3 suggests that the fraction of rotating objects 
does not evolve within this redshift interval. Epinat et al. (2009) 
derive a fraction of rotating galaxies of 22% in a z ~ 1.2-1.6 
sample, but a direct comparison with this value is more difficult 
considering that their definition of rotating objects is quite dif- 
ferent from the one adopted here. 

In any case, these results should be taken with great caution 
since several observational effects may easily affect the capabil- 
ity of identifying rotating objects. More specifically, the lack of 
angular resolution may prevent the detection of a rotation pat- 
tern or, vice versa, may blend two merging systems mimicking 
a rotation curve. Moreover, the limited sensitivity to low surface 



brightness disks (which is a strong function of redshift) may pre- 
vent us to identify large rotating disks in distant galaxies. 

In the following we compare the ordered rotation (the ve- 
locity gradient Av over the map) with the random motions (the 
mean velocity dispersion). 

By using the results of our dynamical modeling we obtain 
the maximum rotational velocity iV,„a.x) and the intrinsic gas ve- 
locity dispersion cr,„, (i.e. the turbulent gas motions). In Fig. 20 
we plot Vmax versus o-jnf for our rotating objects, o-jnf is a con- 
stant value added in quadrature to the model velocity disper- 
sion and determined, as in Cresci et al. (2009), by fitting the 
observed velocity dispersion map. The line Vmax = o-im discrim- 
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Fig. 19. Fit results for the object Q0302-C131. Panels, axis, symbols and color bars indication same as in Fig. 9. 
Table 3. Fit Results. 



Object 



FWHMpsFi"] 



Model parameter 



on 



i n ' 



[/ogio(M,,v„/M0)] 



rAkpcV^' V^aAkms-'] 



SSA22a-M38 


0.5 


149 ±2 


83!? 


SSA22a-C16 


0.5 


-111±2 


60XiS 


CDFS-2528 


0.7 


-25 ±6 


60 (< 70) 


SSA22a-D17 


0.65 


-45 ±5 


80 (> 70) 


CDFA-C9 


0.5 


18±2 


60 (> 20) 


CDFS-9313 


0.7 


128 ±4 


60<5' 


CDFS-9340 


0.7 


44 ±12 


70(5) 


3C324-C3 


0.7 


32 ±7 


60+ '0 


CDFS-14411 


0.6 


78 ±5 


85!^ 


CDFS-16767 


0.55 


98 ±12 


15<5' 


Q0302-C131 


0.3 


48 ±5 


30 (< 60) 



11.34 ± 0.02 (±0.12)<3> 

10.31 ± 0.02 C;5;'J^)(3) 

10.71 ± 0.05 (> 10.65)<'*> 

10.70 ± 0.06 df^P 

9.93 ± 0.03 (!«;^3)<3' 

9.29 ± 0.04 <5' 

10.3 ± 0.1 <^' 

10.12 ± 0.12 c;;*8)(3) 

9.59 ± 0.04 (±0.21)<3' 

10.93 ± 0.11 <^' 
9.99 ± 0.07 (> 9.59)'-" 



2.574 
1.730 
1.644 
1.214 
0.821 
0.713 
1.420 
2.5 
1.240 
0.362 
0.466 



346!^^ 
199+95 

> 197<''> 

260!^^ 
129!2f 

67<5) 
151(5> 

69!!? 
623<5' 

> inc' 



(1) Parameter hold fixed, confidence interval estimated as explained in the text. 

(2) Parameter hold fixed. 

(3) Systematic error. Obtained combining the systematic error contributes due to the inclination estimate and to the PSF width estimate 

(4) Upper-lower limit. 

(5) Parameter unconstrained: we can not give a systematic error contribute due to this parameter. 



inate the loci of "rotation dominated" {Vmax > cri„,) and "disper- 
sion dominated" objects {Vmax < o-im) and, for comparison, in 
local disk galaxies Vmaxlcr ~ 10. From Fig. 20 it is clear that, 
with Vmax/o'int < 2, the majority of the rotating galaxies at z~3 
(6 out of 11) is composed of dynamically "hot" disks, in con- 
trast with local galaxies. Two galaxies have a larger Vmax/o'int 
ratio but still smaller then local galaxies, while three galaxies 
are dynamically "cold". The finding that high-z disk galaxies 
are much more turbulent than local galaxies was akeady ob- 
tained by other studies at z~l-2 (e.g. Genzel et al. 2006, Forster 
Schreiber et al. 2006) and suggestive of high gas fractions mak- 
ing disks gravitationally unstable (Tacconi et al. 2010, Daddi 



et al. 2010a). Forster Schreiber et al. (2009) obtain with the SINS 
sample Vmax/cr ~ 1 - 7 with a mean value of 4.5. Our result at 
z~3.3 {{Vmax/o'i„,)z=3.i ~ 2) goes in the same direction and is 
even more extreme considering that the majority of our galaxies 
have Vmax/o'int comparable with the smallest values observed in 
the SINS survey. We do find however a few galaxies which are 
dynamically "cold" even when compared to local galaxies. The 
missing link between high redshift and local galaxies seems to 
be provided by Epinat et al. (2009) who obtain {Vmax/o'int} ~ 3.5 
(at z ~ 1 .2 - 1 .6) when considering the objects classified as "ro- 
tating disks" and "perturbed rotators". This number rises to 7.2 
when considering only the "rotating disks". 
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Fig. 20. Maximum rotational velocity V,,,^., versus intrinsic velocity dis- 
persion cr„„. Objects from the AMAZE and LSD samples are marked 
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Table 4. Stellar masses and Star formation rates. 
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[togio(M,/Mo)] SFR[Mo/)'] 



SSA22a-M38 
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CDFS-2528 
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In Fig. 21 we compare dynamical and stellar masses for the 
rotating objects. We estimate dynamical masses in the range ~ 
1.8 X IO'^Mq -2.2 X 10"Mo with a mean value of 1.6 x lO'^M©. 
The stellar masses are reported in table 4 and are estimated from 
standard broad-band SED fitting (see Sect. 6 of Maiolino et al. 
(2008a) and Sect. 4.3 of Mannucci et al. (2009) for more details). 

All of the objects are consistent with having stellar masses 
lower than dynamical masses within the \cr uncertainties. By 
ascribing the difference between dynamical and stellar mass to 
the gas mass, the locii of constant gas fraction are given by the 
dashed lines. Here we are assuming negligible contribution of 
Dark Matter, justified by the fact that the spatial region sam- 
pled in our dynamical study is of the order of only 2-3 charac- 
teristic radii of the exponential disk, where the contribution of 



dark matter is still small. In any case, if this assumption is not 
appropriate, we can call this quantity a "gas and dark matter" 
fraction. Three objects require at high confidence a gas fraction 
of 90% or higher. One object (SSA22a-C16) has Httle room for 
any gas mass. Most of the other galaxies have the gas fraction 
poorly constrained or totally unconstrained, except for SSA22a- 
M38 (with a gas fraction constrained within 30% and 80%) 
and CDFS-14411 (with a gas fraction constrained to <40%). 
However, due to the large uncertainties, it is not possible to draw 
firm conclusions on the observed wide range of gas fractions, 
and more data are needed to set tighter constraints on this prop- 
erty. 
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Fig. 21. Comparison of the dynamical and stellar mass for the objects 
analyzed in this paper. The masked gray region represent the non physi- 
cal condition of M^, > M^,.,, . The brown dashed lines represents the loci 
of constant gas fraction (90%, 60% and 30% gas fraction, from left to 
right, respectively). 



It is interesting to compare the amount of gas inferred from 
the comparison between dynamical and stellar mass with the 
amount of gas that can be inferred from the Schmidt-Kennicut 
law. More specifically, if we assume the validity of the Schmidt- 
Kennicut law (hereafter SK) for our z ~ 3 galaxies, we can cal- 
culate the gas surface density from the star formation rate surface 
density, and therefore derive the gas mass. In particular we use 
the relation presented in Mannucci et al. (2009) (eqs. 1 and 2) 
that adopt for the slope of the SK relation the value 1 .4 given 
by Kennicutt (1998). To infer the SFR surface density we use 
the characteristic radius of the mass distribution r^ estimated in 
our dynamical modeling as the radius of the galaxy. Then, under 
the assumption of a negligible contribution in mass of the dark 
matter, and by adding the stellar mass to the gas mass we obtain 
the baryonic mass of the galaxy (M/,,,,). In this way we obtain 
for our objects an estimate of the mass by using only the results 
ofthe SED fitting. 

In Fig. 22 we show the comparison of the Mj,„ estimated 
from the gas dynamics modeling presented in this paper and the 
Mhar estimated from the SK law, combined with the stellar mass. 
The two estimates of the mass, considering the error bars, are 
generally quite consistent. The mean difference (in module) (i.e. 
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the mean scatter of the points around the M^yn - Mha,- condition) 
is of ~ 0.5 dex and there is no particular bias of one of the two 
values on respect of the other The mean residual of the points 
on respect of the Mjy„ = Mhar condition (i.e. the difference in 
module divided by the errors) is ~ 1.2. We can conclude that 
the results of the SED fitting give a reliable estimate of the dy- 
namical mass. This also implies that the two assumptions of the 
validity of the SK law and of a negligible contribution in mass 
of the dark matter appear to be valid. 

Daddi et al. (2010b) using CO observations of galaxy sam- 
ples at low and high redshifts (z ~ 0.5 and z ~ 1.5) find a 
1 .42 slope for the SK relation, very similar to that by Kennicutt 
(1998). Genzel et al. (2010) using systematic data sets of CO 
molecular line emission in z ~ 1-3 normal star-forming galax- 
ies find instead an index of ~ 1.1 - 1.2. In the present work we 
have used the Kennicutt (1998) index for comparison with other 
data in the literature, but using the smaller Genzel et al. (2010) 
index we would have obtained, on average, gas masses larger by 
only ~ 0.3 dex. 
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Fig. 22. Comparison of the dynamical mass M^y,, obtained from the gas 
dynamical modeling and the baryonic mass M^a, obtained by combin- 
ing the stellar mass and the gas mass inferred from the SK law. The 
dashed line represents the condition of equality of the two masses. The 
dot-dashed lines represents the mean scatter of the points around the 
Mrfvn = M],a,- condition (~ 0.5 dex). 



5.2. The Tully-Fisher relation at z ~ 3 

The Tully-Fisher relation (hereafter TFR) is the relation between 
luminosity (or stellar mass) and maximum rotational velocity 
of disk galaxies (Tully & Fisher 1977). This relation, originally 
used as a distance indicator for disk galaxies, represents an im- 
portant feature for understanding galaxy formation and evolution 
because it links directly the angular momentum of the dark mat- 
ter halo with the stellar mass content of disk galaxies. In fact, 
according to theoretical models, galaxy disks form out of gas 
cooling down from a surrounding dark matter halo, maintain- 
ing its specific angular momentum and settling in a rotationally 
supported disk. Therefore, the structure and dynamics of disk 



galaxies at different cosmic epochs are expected to be closely 
correlated with the properties of the dark matter halo in which 
they are embedded. 

In recent years, an increasing number of dynamical studies 
of disk galaxies in the local and intermediate-redshift Universe 
allowed the investigation of the evolution of the TFR with red- 
shift. 

The evolution of the TFR is expected to be related both to 
the conversion of gas into stars and to the inside-out growth 
of the dark matter halo by accretion. In fact, while the extent 
of the halo grows considerably with time, the circular veloc- 
ity of the halo grows less, keeping the rotation curve approxi- 
mately flat to larger and larger distances. The accretion of the 
dark matter is followed by accretion of baryonic gas, which is 
subsequently converted into stars by ordinary star formation in 
the disk. However, the details of the process and the amount of 
evolution expected depends strongly on the model assumptions, 
on the accretion mechanism and the timescale needed to con- 
vert the gas into stars. Any successful model of disk formation 
should then be able to reproduce the slope, zero-point, scatter, 
and redshift evolution of the TFR. 

There are still discrepant results on a possible evolution at in- 
termediate redshifts of the tight relation observed in local galax- 
ies (e.g., Haynes et al. 1999; Pizagno et al. 2007). For example 
Vogt et al. (1996) reported very little evolution of the B-band 
TFR up to z ~ 1, but other groups (e.g. Harden et al. 2003 and 
Nakamura et al. 2006) found a strong brightening of ~ 1-2 mag 
in B-band luminosity over the same redshift range. 

Furthermore, the interpretation of a possible evolution of the 
luminosity TFR is difficult since the luminosity and the angu- 
lar momentum might be evolving together. As a consequence, 
the stellar mass TFR, which correlates the stellar mass and the 
maximum rotational velocities of disks, offers a physically more 
robust comparison as it involves more fundamental quantities. 

In this context, Flores et al. (2006), Kassin et al. (2007), and 
Conselice et al. (2005) have found no evolution in both the slope 
and zero point of the stellar mass TFR up to z = 1.2. 

In contrast Puech et al. (2008) detect an evolution of the TFR 
zero point of 0.36 dex between z ~ 0.6 and z = 0, by using a 
sample of 18 disk like galaxies observed with the integral field 
spectrograph GIRAFFE at the VLT. 

At z~2 Cresci et al. (2009) build a TFR at z ~ 2 based on 
the SINS sample and find an evolution of the TFR zero point of 
~ 0.4 dex relative to the local TFR. 

By using the dynamics of z ~ 3 galaxies inferred from our 
data, we are now able to derive the TFR at z ~ 3. In Fig. 23 we 
present the z ~ 3 TFR obtained by using the data presented in 
this paper (obviously only for the rotating objects). The maxi- 
mum rotational velocities have been computed from the best fit 
models as reported in Tab. 3. The stellar masses are those re- 
ported in table 4. 

All our data points are consistent with having a stellar-to- 
dynamical mass ratio smaller than the local value. We can actu- 
ally observe an evolution of the TFR with the redshift and this 
evolution is in the direction predicted by the models: the dynam- 
ical mass is already in place at this cosmic epoch, but the stellar 
mass has yet to be formed. 

The large scatter (~1.5 dex) of the TFR at z ~ 3 suggests 
that at this redshift the relation is not yet in place, probably due 
to the young age of the galaxies at this epoch (between 50Myr 
and SQQMyr). However we can not exclude that our V„,ax mea- 
surements are affected by an incorrect assumption of a rotating 
disk kinematics for some objects of our samples, which intro- 
duces a large scatter in the relation. 
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Fig. 23. Tiie Tully-Fisiier relation at z ~ 3 reconstructed by using the 
dynamical and stellar masses of the rotating galaxies found in our sam- 
ple. The solid black line represents the local relation, the magenta line 
represents the z ~ 2 relation (Cresci et al. 2009) and blue solid line rep- 
resents the fit to our data (z ~ 3) by keeping the same slope as for the 
local and z ~ 2 cases. The dotted blue lines represents the formal Icr 
uncertainty on the zero point of our fitted relation. 
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Fig. 24. The 5 0.5 TuUy-Fisher relation at z ~ 3. The solid black line rep- 
resents the z ~ 1 relation (Kassin et al. 2007), the dashed line represents 
the fit to our data (z ~ 3) by keeping the same slope as for the z ~ 1 
relation. 



Assuming that the TFR is in place, despite the large data 
scatter, and that the slope is the same of the local universe, we 
can fit the relation obtaining a zero point for our data of -0.97 + 
0.13. This implies an evolution of the zero point of 1.29 dex 
relative to the local TFR and 0.88 dex relative to the z~2 TFR. 



Given the high intrinsic velocity dispersion observed in our 
z ~ 3 sample, we also explored the 5 0.5 - VO.SV^ + cr^ esti- 
mator introduced by Kassin et al. (2007) to study the impact of 
turbulent motions. In Fig. 24 we present the 5 0.5 TFR obtained 
by using our Y,nax and cr,„, values. Kassin et al. (2007) found a 
tight relation (~ 0.3 dex scatter in M*) in a sample of galax- 
ies with redshift up to z = 1.2, with no significant evolution. 
By only allowing for a change in zero point of the z ~ \ rela- 
tion we estimate a mean scatter of ~ 0.5 dex in our galaxies at 
z ~ 3.3. The reduction of the scatter in the S0.5 TFR compared 
to the classical one (~ 1.5 dex) suggests that, at least for some 
objects, part of the mass might be supported by turbulent or not 
ordered motions. We observe a significant evolution of the zero 
point of the relation between z ~ 1 and z ~ 3 (~ 0.8 + 0.1 dex 
in M*), providing further support for the observed evolution on 
the classical TFR. 

The use of SINFONI integral field spectra represents a major 
step forward in the modeling of the dynamics of high-z galax- 
ies, as the two-dimensional mapping is more complete and not 
biased by a priori assumptions about the kinematic major axis 
and inclination of the system as in long-slit spectroscopy. Also 
SINFONI data can allow a better selection of the rotating sys- 
tems and ensure a much lower contamination from complex 
dynamics and irregular motions than in longslit data sets. This 
might be a reason for the apparent inconsistency of the results of 
kinematical studies at redshift greater than 0.5. 

6. Conclusions 

We conducted an extensive investigation of the dynamics of 
z ~ 3 galaxies using near infrared integral field spectra obtained 
with SINFONI at the ESO VLT. This is the first time that such 
dynamical analysis is conducted for a relatively large sample (33 
galaxies) at z ~ 3. 

- We estimated galaxy masses by fitting rotating disk mod- 
els to selected galaxies. Due to the relatively low S/N of the 
data, we implemented a simple but quantitative method to 
identify smooth velocity gradients, which are the signature 
of possible disk rotation. After this pre-selection, we per- 
formed complete dynamical modeling using a rotating disk 
model with an exponential mass profile. 

- We found that the kinematics of ~ 30% of the objects (11 out 
of 33) is consistent with that of a rotating disk. This fraction 
is fully consistent with that found with z ~ 2 samples, like 
SINS, suggesting that the fraction of rotating objects does 
not evolve between z ~ 2 and 3. 

- When considering only the galaxies observed with AO and 
finer pixel scale (LSD sample), only 1 out of 9 shows a 
smooth velocity gradient. This is likely the consequence of 
such observational setup providing higher spatial resolution 
but lower sensitivity to the outer, low surface brightness re- 
gions which dominate the "rotation signal". 

- The 1 1 rotating objects have dynamical masses in the range 
~ 1.8 X 1O''M0 - 2.2 X 10"Mo with a mean value of ~ 
L6xlO'°Mo. 

- By comparing the amplitude of the inclination-corrected ve- 
locity gradient with the intrinsic gas velocity dispersion de- 
rived from models, we found that the majority of the rotating 
galaxies at z ~ 3 are dynamically "hot" disks with an aver- 
age {Avdeproj/cTini) ~ 0.9. This Confirms the results obtained 
by other studies at z ~ 1-2 (Av/cr ~ 2 - 4), indicating that 
high-z disk galaxies are more turbulent than local galaxies 
and even more "hot" than galaxies at z ~ 2 . 
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- We compared stellar and dynamical masses for the rotating 
objects, and derived an estimate of the "gas and dark mat- 
ter" fraction. Using independent estimates for stellar mass 
and star formation rate, and applying the Schmidt-Kennicut 
law, we obtained a SED-based estimate of the total baryonic 
mass of our objects. This estimate is consistent with our esti- 
mate of the mass from gas dynamics (rms scatter lower than 
O.Sdex). This shows that SED fitting gives a reliable estimate 
of the dynamical mass and also implies that our assumptions 
on the validity of the SK law and of a negligible contribution 
of dark matter are well justified. 

- Finally, we obtained the Tully-Fisher relation at z ~ 3. All 
our data points are consistent with a stellar-to-dynamical 
mass ratio smaller than the value in the local universe con- 
firming the redshift evolution of the relation already found at 
Z ~ 2. This is consistent with models predictions, according 
to which stellar mass is still building up compared to the to- 
tal dynamical mass. However, the large scatter of the points 
may also suggest that, at this redshift, the relation is not yet 
in place, probably due to the young age of the galaxies. The 
Tully-Fisher relation based on the 5 0.5 indicator has a signif- 
icantly smaller scatter that the "classical" one and also shows 
an evolution in zero point compared to z ~ 1 . The reduction 
in scatter suggests that, at least for some objects, part of the 
mass might be supported by turbulent motions. 
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